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power density can be reduced by a factor of four as com- 
pared to a single horn feed. Additional reduction is possible 
by proper selection of feed defocusing. 
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On Thin-Wire Multiturn Loop Antennas 


CLAYBORNE D. TAYLOR, MEMBER, IEEE, AnD CHARLES W. HARRISON, JR., reLLow, IEEE 


Abstract—-A general theory of thin-wire multiturn loop antennas 
is presented. For convenience the windings of the loop are considered 
to form a circular helix. An integral equation is derived and solved 
numerically for the current distribution on the antenna. The antenna 
impedance, efficiency, and pattern are obtained. 


INTRODUCTION 


HE MULTITURN loop antenna has important ap- 
plications as both a transmitting and a receiving 
antenna, primarily because the volume occupied by the 
antenna is not required to have dimensions comparable 
to the operating wavelength. However, the multiturn loop 
is formidable to analyze because of its somewhat complex 
physical configuration. In particular it is difficult to deter- 
mine the high frequency behavior. In this paper a multi- 
turn loop in the form of a circular helix is considered, and 
a general theoretical-numerical formulation is presented. 
The single-turn loop antenna has been studied by a 
number of investigators [1 ]-[6]. Generally a delta function 
generator is used for the driving mechanism in the analy- 
sis; however, in the formulation presented here a magnetic 
ring source is used. Recently the polygon loop antenna, 
single turn and multiturn, was studied by using the so- 
called reaction technique developed by Richmond [7 ], [8 ]. 
By using certain symmetry properties and judicious 
approximations, Padhi [9] has developed a theory of 
multiturn loop antennas with particular winding configu- 
rations. Recently, Smith [10] developed a formulation for 
electrically small multiturn loops. Also, Smith considered 
proximity effects on the circumferential current distribu- 
tion about the wire element cross section. The research 
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reported here is not restricted to electrically small loops 
and is more general than that presented by Padhi. 


ANALYSIS 


Consider the multiturn loop antenna to be in the form 
of a circular helix with the following parametric equations: 


x = Ry cos é 
y = Rysiné 
(1) 


where @ is the parameter, Ry is the radius of the helix, and 
p is the pitch of the helix. The parameter @ may be ex- 
pressed in terms of arc length along the helix; ie., 


6 = as (2) 


z= pe/2r 


where 
a = [Re + (p/2r)? 7. (3) 


The vector and scalar potentials, when evaluated on 
the axis of the wire forming the multiturn loop, may be 
written [11] 


Ae af (8-8) I.(s!)G(s,8') (4) 
60) - 52, fw SUC) 6) 

where 
G(s,s’) = exp [—jk(R?(s,8’) + a”)? ]/(BR2(s,8") + a?) ¥? 
(6) 
R?(s,s') = (2Ro)? sin? [a(s — 8’) /2] + (ap/2r)*(s — 8’)? 
(7) 
$5 SUR eel = ade Ge. (8) 


In these expressions a is the radius of the wire forming 
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the multiturn loop and 7 = (n/c)? is the intrinsic wave 
impedance of the surrounding medium. 

On the axis of the wire the component of the total elec- 
tric field tangent to the axis is 


E,(s) = E,°(s) — ~ 4(9) —jwoA(s)=0 (9) 


where F,i"¢ is the tangential component of the incident 
electric field (or an impressed electric field as in the case 
of a driven antenna) and the potential functions represent 
the scattered field. Substituting (4) and (5) into (9) 
yields 


) a 
i ds' — I,(s') — G(s,s") + wf ds’ I,(s')§-8'G(s,s') 
Lt 9s as L 


Ark 
=-j ee a (10) 


an integrodifferential equation for the current distribu- 
tion. Observing that (6) and (7) indicate 


a ? a ? 
. G(s,s’) = -— aa G(s,s’) 
and integrating by parts (10) becomes 


fas |4 “ i058) I,(s')G(s,s') 
os 


d | 7 

= gai Be(8') G(8,8") : 

If the wire of the loop is lossy, an internal impedance 
per unit length z* may be defined 


E,™°(s) losses = —2'I(s) (12) 


which is added to the expression for the incident electric 
field component derived in Appendix A for the multiturn 
loop antenna driven symmetrically from both ends by 
magnetic ring sources. The radii of the magnetic rings 
are chosen to represent the actual source driving the loops 
(see Appendix A). 

Inherent in the foregoing analysis are certain restric- 
tions. For the thin-wire approximations to be valid, the 


Ark _. 
= ae (11) 


radius of the wire helix must be much greater than the ~ 


wire radius, and the wire radius must be much less than a 
wavelength. In this analysis the proximity effect on the 
circumferential current distribution is ignored; therefore, 
the inequality p > 4a must be satisfied [9], [10]. 

Evidently the wire connecting the ends of the helix does 
not conform to the helical geometry. If the separation of 
the ends is small in terms of wavelength, a very simple 
procedure for obtaining the contribution of currents and 
charges on the wire connecting the ends may be used (see 
Appendix B). 


NUMERICAL SOLUTION 


A convenient numerical solution of the integrodifferen- 
tial equation for the current distribution may be effected 
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by the so-called method of moments [4]. A piecewise 
sinusoidal function is used to represent the current dis- 
tribution as opposed to using the, perhaps more common, 
flat pulses or triangle pulses [7]. The current distribution 
is then 


N 
I,(s’) = > In(s’) U (8' j8my1)Sm) (14) 
m=1 
where 
fate). CE See) eee) oa) 
sm K(Smn41 es Sm) 
1, Sm < s’ < Sm41 


U(8' s8m41,S8m) = 
0, otherwise 


and where {s,,} is an ordered set of N + 1 points along 
the axis of the wire including the end points. Further- 
more, it is required that there be no discontinuity in the 
current at the wire ends. 

Taking the derivative of the current distribution yields 


a, La eens a E T ’ Je ,, 1 
a 3(s’) -z om m(8’) | U(8';8m41,8m). (16) 


Operating on the current distribution with the Helmholtz 
operator yields 


& 
— + 2 f 
Fe : | ee 


Nd 
= DS In(s)[6(5! — bn) ~ 88! — sma) 


d 
— a I,(s’)[6(s’) — 6(s’ — L)]. (17) 
By substituting (17) into (11) one obtains 
N d 89S atl 
>, | In(s’) G(s,8’) 
m=1 ds sf=sm 
Smt) 
_ we f ds’ (1 — 8-8’) Im(s’) G(s,s’) 
Aok | 
—j i ZI n(s) U(sitmnm) | 
v7] 
Ark 
= a E,'**(s). (18) 


A linear system of equations is obtained from (18) by 
the introduction of a convenient set of points on which 
(18) is satisfied. Also, the boundary conditions on the 
current must be enforced. For the multiturn loop the 
current boundary condition is 


I,(0) = I,(L) (19) 


TAYLOR AND HARRISON: MULTITURN LOOP ANTENNAS 


or 
(20) 


a = aN41. 
Using (15) and (20) in (18) and satisfying (18) at points 
S = Sq, for n = 1,2,+++,N, one obtains 


N 
> Tamm = Tr, 


m=1 


(21) 


where 


Tha = [G (Sp,8x ) = LG (Sn81) + G(8a,8v41) | 
miei 
eos: Avs Saga) ksin kA, 


52 
_ J ds’ (1 — 8:8’) sin k(s2 — s’)G(8n,8') /sin kAy 
0 


L 
fe / dl = $53) aE = G8 Jan BA 
aN 


Anz 
— J— bu. (22) 
nk 
For m > 1 
G(Sns8m—1) — G(Sn;8m) COS KAn—1 
Tam ae : 
k sin kAn_1 
4 G(Sn,8m41) — G@(Sn,8m) COS KAm 
k sin kA,, 
+f ds’ (1 — 8°38’) 
8m-1 
-sin h(Sm—1 — 8’) G@(8n,8’) /sin kKAy_1 
8m+t 
- ds’ (1 — 8,°8’) 
; ; Anz 
sin K(smi1 — 8’)G(Sn,8’)/sin kA, — j ps Sam: 
n 
(23) 
For n = 1,.N 
4 Fy,ine a 
r= —-j Aa Bie°(Sn) (24) 


nk 


Here An = Smyi ~— Sm, for all m. It may be noted that 
satisfying (18) at s = svi1 = L yields the same equation 
as satisfying it at s = s, = 0, provided that (20) is used. 


APPLICATIONS 


The characteristics of interest concerning the loop an- 
tenna are its impedance, efficiency, and pattern. Since 
the loop is driven from the wire ends with voltage V, the 
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(7,8, 6) 


Fig. 1. Multiturn loop in form of circular helix shown with co- 
ordinate system. 


impedance is (see Appendix A and Fig. 1) 


Za = aes (25) 
AO) TAL). 
The antenna efficiency is 
2] 


~ GRe [Z.]| 1.(0) [7]. 


For a solid or tubular conductor with a wall thickness 
that is large compared to the skin depth the internal 
impedance per unit length is [1] 


2zé= (1+ 7) (2rac6,)—. (27) 
The skin depth 6, is given by 
5s = (2/wpo)¥?, (28) 


To determine the field radiated from the loop, the re- 
striction that the ends of the loop be separated a distance 
small compared to the wavelength is again imposed. Ac- 
cordingly, the radiation pattern is essentially that of a 
single-turn loop with the current 


NT 
T(g') = DLL! + 2(n — 1a} /a]. (29) 


n=1 


Here NT is the number of turns and @ is given in (3). 
Using (29) and standard techniques, the electric field 
components in the radiation zone are obtained in the fol- 
lowing form: 


aV_ exp (—jkr) 


QnDna Po (4,6) 


Ee (7,0,¢) = —j 


nV exp (— jkr) 
r 


QnZa Fs (9,¢) 


E4(7,0,0) = —j 
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where the pattern functions are 


k A 2a 
F.(6,8) = io [ae 108 
-exp [ jkRo sin 6 cos (@ — ¢’) ] sin (@ — ¢’) 
F5(0,6) = ; a ae dg! I(¢!) 


-exp [ jkRo sin 6 cos (@ — ¢’) ] cos (6 — ¢'). 


The impedance of a multiturn loop antenna operating 
in the low frequency regime as given by King [14] is 


Za = Ra + 7X (30) 
where 
Re =F n( ko) (NT)? + rity + C(QeRat) + oP NT 
(31) 
8Ro a 
X, = nkRo {wr jm —- 2| +2 x (NT — ano} 
(32) 
u(t) = — 2 Ble) +l = =e; = — «,) K(e) 
es = [1+ p?/ ee 


with K(e:) and E(e;) as complete elliptic integrals of the 
first and second kinds, respectively, and r; the internal 
resistance per unit length for the wire forming the loop 
antenna. The foregoing formula attributed to King has 
been simplified somewhat and extended to include losses. 
Further, it is modified in accordance with the geometric 
configuration considered here being different from that 
considered by King. 


NUMERICAL RESULTS 


First the impedances of 1-, 2- and 3-turn perfectly con- 
ducting loop antennas are shown in Table I. The im- 
pedances are compared with those from the low frequency 
approximation and the results from King and Harrison 
[6] for single turn loops. In general the agreement is 
quite favorable. At low frequency the real part of the 
integral equation solution is not as accurate as the low 
frequency approximation because of accumulative trun- 
cation errors in the numerical computations (the computer 
used to obtain the presented data is the UNIVAC 1106 
that carries 8 significant figures in single precision com- 
putations). For example if the reactance is accurate to 
+0.1 @ then the resistance is probably accurate only to 
within 0.1 @ In any event the foregoing comparisons 
are sufficiently good to indicate the validity of the assumed 
driving mechanism and the computer program. 

From the foregoing data it is seen that adding turns to 
the loop antenna not only lowers the resonant frequency 
but it also causes the resonances to become much sharper. 
This characteristic is also exhibited in Figs. 2 and 3, 
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TABLE I 
Mou.trrurn Loop Imeepance 
KR, 2a 
NT = 1 NT = 2 NT = 3 
0.05 0.018 + j 60.0 0.036 + 4 166.3 0.057 + 4298.9 


(0.001 + j 61.1)* (0,005 + j 193.6)* (0.0111 + 5 374.5* 


0.10 0.087 + j 124 0.323 + j 477.1 8.75 + 43502 
(0.020 + j 122)* (0.079 + j 387.2)* 
0.15 0.272 + 4 197 18.1 + § 3112 0.35 - 4 627.6 
(0.10 + j 283)* 
(0.139 + j 203)** 
0.20 0.762 + 5 287 2.58 - j 989.0 0.082 ~ 4 273.8 
(0.316 + 5 244)* 
(0.535 + j 296)** 
0.25 2.38 + 4 454 0.718 - j 432.6 0.026 - j 141.1 
(1.82 + j 422)** 
0.30 6.84 + 4 656 0.386 - 4 260.1 0.019 - 4 53.0 
(6.11 + 4 614)** 
0.35 23,1 + 4 1018 0.241 - j 166.4 0.107 + j 29.0 
(23.2 + 3 967)** iu 
0.40 118 + 4 1929 0.147 - j 101.1 0.637 + 4 136.6 
(134 + 5 1901)** 
0.45 4023 + 4 8492 0.085 - 4 47.9 5.19 + 5 388.1 
(7942 + 4 9076)** 
Q = 2 &n(20R fa) = 10.0, p = 0.0258 m 
Ro = 0.1524 m, zi = 0 
*obtained from low frequency approximation 
*kobtained from King and Harrison [6_] 
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Fig. 2. Input resistance of a 5-turn copper she loop antenna. 
Ry = 0.2 m, a = 0.00079 m, and p = 0.01 m. 
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Fig. 3. Input reactance of 5-turn copper wire loop antenna. Ry = 0.2 


m, p = 0.01 m, and a = 0.000794 m. 
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where the antenna resistances and reactances are shown 
for a 5-turn copper wire loop antenna. The experimental 
results were obtained by Flaig [15]. For the calculated 
results the number of piecewise sinusoidal sections used 
is 82 while utilizing the symmetry of the antenna current 
distribution. The convergence of the numerical solution 
depended upon the wire radius as reported previously [18 ]. 
For that reason the computer program could not handle 
as thin a wire as Flaig used to construct his 5-turn loop 
antenna. The parameters of Flaig’s antenna are Ay = 
0.2 m, a = 0.000794 m, and p = 0.01 m; the parameters 
used in the computer program are R = 0.2m, a = 0.005m, 
and p = 0.02 m. In both cases the same internal impedance 
per unit length of copper wire is used. The agreement 
between measured and calculated input reactances is 
quite good. However, the input resistances do not com- 
pare nearly as well; yet the agreement tends to improve 
at higher frequencies. 

In making measurements Flaig considered his balum 
transformer (Anzac Model XT-617) behaved essentially 
as a transmission line (with phase velocity of free space) 
in series with a 4:1 impedance transformer. How these 
assumptions affected the antenna impedance measure- 
ments is not known. However, at low frequencies the 
experimentally determined transmission line characteristic 
impedance departed significantly from the assumed value 
based upon an average measured value for the frequency 
range from 50 to 150 MHz. The first resonance occurs 
well below 50 MHz. Furthermore, the low frequency ap- 
proximation agrees better with the numerically obtained 
results. For example at f = 12.5 MHz the low frequency 
approximation yields Z, = 1.34 + 71539; the integral 
equation solution is Z, = 1.52 + 71059. Thus it is ex- 
pected that experimental error accounts for a large part 
of the disagreement between the measured and calculated 
input resistances. 

Antenna efficiency is an important consideration for 
the multiturn loop. Flaig used an analysis developed by 
Munk [16] to compute the efficiency as well as the radi- 
ation resistance (Munk did not obtain the reactance). 
Theoretical and experimental results are shown in Fig. 4 
for a 5-turn copper wire loop. At low frequencies, even 
past the first resonance, the antenna efficiency of the 
5-turn loop is poor. But near the second resonance the 
efficiency tends to improve. 

Field patterns of the 1-, 2-, and 3-turn loop antennas 
are shown in Figs. 5 and 6. The patterns are for conditions 
near the first resonance of the antennas. Apparently the 
patterns of multiturn loops below the first resonance are 
similar to the low frequency pattern of a single turn loop. 


APPENDIX A 


For convenience the loop antennas are considered to be 
driven from magnetic current sources. These magnetic 
sources may represent the equivalent field distribution in 
the driving gap region such as the open end of a trans- 
mission line [17]. Recently Tsai [12] has derived expres- 
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Fig. 4. Antenna efficiency of 5-turn copper wire loop. Ry = 0.2 m, 
a = 0.000794 m, and p = 0.01 m. 
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Fig. 5. Field patterns of 1-turn, 2-turn, and 3-turn loop antennas. 
Q = 21n (27R)/a) = 10.0, p = 4a, and N = 50. (For these data 
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Fig. 6. Field patterns of 1-turn, 2-turn, and 3-turn loop antennas. 
Q = 2)n (27R,/a) = 10.0, p = 4a, and N = 50. (For these data 
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sions for the field about two magnetic current sources: 
the magnetic frill and the magnetic current ring. The 
magnetic frill is the magnetic current equivalent to the 
TEM mode field distribution in the gap of an open-ended 
coaxial transmission line connected to a perfectly con- 
ducting plane and the magnetic ring is the magnetic cur- 
rent equivalent to the field distribution in an infinitesi- 
mally thin circumferential slot driving a cylindrical an- 
tenna. 

One of the advantages of using the aforementioned 
magnetic current sources is that analytic expressions are 
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available for the field about the sources. For example the 
field along the axis of a magnetic ring is 


a Va a K (2,20) 


E.(p,2,0) | po _ 2 da 


(Al) 


where a is the radius of the ring and 


exp [—jk((z — 29)? + a?) ¥] 
((z — 20)? +- a?)2 


with z locating the center of the ring along the z axis and 
V the voltage impressed across the equivalent infinitesimal 
gap. It is straightforward to extend (A1) to represent the 
axial electric field produced by the equivalent magnetic 
currents for a general gap excitated antenna. The result is 


K(z,%) = 


E.(,2,6) | 0 = + | 


E K(e2) | Ee? (z') dz’. (A2) 
gap LOG 

However, the foregoing result requires a determination of 
the electric field impressed in the gap. Unfortunately this 
electric field is not generally known a przorz. 

When (A1) or (A2) is used for the impressed field the 
integral equation for the current distribution is obtained 
by requiring the total z-component of the electric field to 
be zero along the wire axis. This procedure is the applica- 
tion of the so called extended boundary condition as dis- 
cussed by Taylor and Wilton [18]. The use of (Al) and 
the extended boundary condition should be mathemat- 
ically equivalent to the use of the delta gap feed and the 
boundary condition on the tangential component of the 
electric field at the antenna surface. The latter procedure 
encounters difficulties when a colocation solution technique 
is used such as is presented in the foregoing. 

At large distances from the gap, ie., 2 >> gap width, 
(A2) yields that the electric field produced by the 
equivalent sources for a finite gap is essentially the 
same as that produced by an infinitesimally thin mag- 
netic current ring. However, at points near the gap where 
k( (2 — 2)? + a@)¥2 < 0.1 


3 a 1 
E.(p,2,9) | o=0 = = a 
gap 


ke Ul , 
aes al E(2') dz 


i 2 
+5 | E,s°?(2') de’ (A3) 


gap 
where R = ((z — 2’)? + a®). A convenient. approxima- 
tion for the electric field in the gap is! 


E f°? (z) ~— = 
g 


where w, is the gap width, so that (A3) becomes for 
a=0,2=0 

V 1 

2 (a + (u_/2)). 


(A4) 


Li, (p,2,0) | p=0; z=0 = 


1 This approximation is well known to yield good antenna im- 
pedances [19]. 
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For a infinitesimally thin magnetic current ring with radius 
a; the axial component of the electric field at the center of 
the ring is V/2a,. Thus (A4) may be used to define an 
equivalent radius of a infinitesimally thin ring to repre- 
sent a band of magnetic current. This equivalent radius is 


as = (a? + (w/2)?)¥?, 


The multiturn loop is then considered to be driven from 
both ends by equivalent magnetic current rings. This is 
a reasonable alternative to driving the loop from the 
center of the wire connecting the ends of the windings 
since the separation of the ends of the windings is an 
electrically short distance. Unfortunately when consider- 
ing the loop antenna, knowing the field along the axis of 
the magnetic ring is not sufficient. In this case the electric 
field component tangent to the loop axis must be known. 
There are two available routes to be taken. First, one 
may numerically compute this field component as Tsai 
did, or one may judiciously obtain an approximation that 
is convenient for computations. The latter approach will 
be followed. 

At distances that are large in comparison with the radius 
of the magnetic ring, the field distribution is essentially 
that of an electric dipole with moment 


P = ra, 2eVi. 
The electric field about a dipole is 
= Pe) (2 


7 


E(r) +34) [30(7-P) — PY 


Ate 


— "tx @x PI] (As) 


where r is the radius vector from the dipole to the point 
of observation. (It is readily noted that (AS) has the 
same far field behavior as obtained by Tsai, provided 
ka <1.) Furthermore, it is observed that (A5) reduces to 
(Al) along the z axis when r—> R, = (7? + a,2)"2, Even 
with r—> R,, (A5) still exhibits the proper far-field be- 
havior. Therefore, the component of the illuminating 
electric-field tangent to the wire loop at point s is approx- 


imately 
eH) yy (2 4 j ® 
gE ae tae 


-[3 (78) (#2) — 8-2] — = LF x (7 X 2] 


E,i?*(s) = 


(A8) 


where § is the unit vector tangent to the axis of the wire 
forming the loop and R, = (R*(s,%) + a2)? with R 
given by (7). Note that (A6) is precise very near the 
magnetic ring R(s,s) ~ a and at distances R(s,s9) >> a’. 
For a single turn circular loop driven from a magnetic 
ring source at s) = 0, the incident field component driving 


TAYLOR AND HARRISON: MULTITURN LOOP ANTENNAS 


the wire loop is 


a2V 1 : 1 as 
ine ee ees . —i{1 ar 
ee ar Fes E e zl as sl 
Ke. as . 
~ oR, =| SENS eH AE) 
where 


Rs = ((2Ro sin as/2)? + a2)". 


Using (A7), results are obtained that agree to within a 
few percent with those presented by King and Harri- 
son [6]. 


APPENDIX B 


In this appendix a procedure is developed to account 
for the wire segment connecting the ends of the helical 
loop. The contribution to the incident field from the cur- 
rent and charge sources on the subject wire segment may 
be expressed. 


#£,"(s) = 


AL 
. 7 0 0 
ees dz’ I, , Y ; t 
ere Zz oe (2') re G,(s,2’) 


+k? (8-2) L.(2')Gi(s,2’) (BI) 


where the wire segment is parallel to the z axis and extends 


from 2’ = 0 to 2’ = AL = pNT. The number of loop 
turns is designated by V7. Also, 


(B2) 


and 
G1 (8,2) 


exp L[—jk((2Ro sin as/2)* + (aps/2mr — 2’)? + a®)V2] 
((2Ry sin as/2)? + (aps/2r — 2’)? + a?)i2 


(B3) 
After some mathematical manipulation (B1) becomes 
4 AL Fe) 2Z=AL 
Ee =j—|]— — I1.(2’)G\(s8,2’ 
(3) J mal L 02’ “ ) ls 4 ) 2f/=0 
AL f°" of @ c ‘ : 
ae [ dz (3 + i) L,(z )| Gi(s,2') 


: ne re 10 
— ak? sin os f df —higji= 
0 0z 


ROR 
[2 (—jhR) ] 
R 
where 


FR? = [2Ry sin (as/2) ? + (aps/2m — 2’)? + a’. 


(B4) 


(B5) 


Over the short-wire segment the current distribution 
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may be considered sinusoidal; i-e., 
I,(Q) sin k(AL — z) + I,(AL) sin kz 


z = ; 6 
T(z) sin KAL cB6) 
From the symmetry of the excitation 
I,(0) = I,(AL) = I,(0) 
then ; 
in k(AL — in kz 
I.(2) = L.(0) sin k( z) + sin (B7) 


sin kAL 


Substituting (B7) into (B4) results in the vanishing of 
the first integral of (B4). The resulting form for /£,°" is 
then added to the incident field in the integral equation 
for the antenna current distribution. Without this correc- 
tion, meaningless results are generally obtained for the 
antenna current distribution. 
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